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►> \ Abstract 

0\ ■ We study some structural aspects of the evolution equations with Pomeron loops recently 

derived in QCD at high energy and for a large number of colors, with the purpose of clarifying 
their probabilistic interpretation. We show that, in spite of their appealing dipolar structure and 
of the self-duality of the underlying Hamiltonian, these equations cannot be given a meaningful 
interpretation in terms of a system of dipoles which evolves through dissociation (one dipole 



splitting into two) and recombination (two dipoles merging into one). The problem comes from 
the saturation effects, which cannot be described as dipole recombination, not even effectively. 
We establish this by showing that a (probabilistically meaningful) dipolar evolution in either the 



target or the projectile wavefunction cannot reproduce the actual evolution equations in QCD. 
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1 Introduction 



Over the last year, an intense activity in the field of high-energy QCD has been triggered 
by the observations that (i) the gluon number fluctuations at low energy play an im- 
portant role in the evolution towards gluon saturation with increasing energy [1-3], and 
(ii) the relevant fluctuations are in fact missed [4] by the existing approaches to non- 
linear evolution in QCD at high energy, namely, the Balitsky-JIMWLK equations [5-8]. 
To cope with that, new equations have been proposed [4,9-11] (see also Refs. [12,13]), 
which include the particle number fluctuations together with the saturation effects in the 
limit where the number of colors iV c is large. At the same time, an ambitious program was 
launched, which aims at generalizing these equations to arbitrary values of N c [14-21]. 
This effort led already to some important results — in particular, the recognition [11, 15] 
of a powerful self-duality property of the high-energy evolution, and the construction of 
an effective Hamiltonian which is explicitly self-dual [17, 20] — , but the general problem 
is still under study, and the evolution equations for arbitrary N c are not yet known. 

In this paper, we shall restrict ourselves to the large-iV c limit and investigate whether the 
evolution equations proposed in Refs. [4, 9] can be given a probabilistic interpretation in 
the sense of the dipole picture. To motivate our search, let us remind here that most of 
the modern approaches to high-energy evolution in QCD have an underlying probabilistic 
structure, which makes their physical interpretation transparent and greatly simplifies the 
analysis of their consequences, via analytic or numerical techniques. 

Specifically, the color dipole picture introduced by Mueller [22-24] describes the BFKL 
evolution [25] at large N c (including the particle number fluctuations alluded to above) as 
a classical Markovian process in which a system of dipoles evolves through dipole splitting. 
This manifest probabilistic interpretation has enabled Salam to build up a Monte-Carlo 
code allowing for numerical studies of the dipole evolution [26]. The dipole picture is also 
underlying the derivation by Kovchegov [27] of the non-linear Balitsky-Kovchegov (BK) 
equation. Recently, this derivation has been extended, by Levin and Lublinsky [28], to 
the large-7V c version of the Balitsky-JIMWLK hierarchy. 

Similarly, the color glass condensate formalism [7,29,30], which aims at a description of 
the high density environment at, or near, saturation, is explicitly formulated in terms 
of a classical probability distribution, which evolves with increasing energy according 
to a functional Fokker-Planck equation, the JIMWLK equation [6-8]. The basic micro- 
scopic process at work here (in addition to BFKL evolution) is gluon recombination : the 
JIMWLK Hamiltonian includes vertices for n — > 2 gluon merging with n > 2, which lead 
to gluon saturation at high energy. The stochastic nature of the JIMWLK evolution is 
most clearly exhibited by its reformulation as a Langevin equation describing a random 
walk in the space of Wilson lines [31]. This is also the most convenient formulation for 
numerical studies, as demonstrated by the analysis in Ref. [32]. 

Note that the color dipole picture can be also formulated as a "color glass" [10, 24], but the 
ensuing picture is complementary to the JIMWK evolution: it describes gluon splitting (at 
large- N c ), but not also gluon recombination. Clearly, a complete picture should involve 
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a synthesis of these two approaches, but such a synthesis has not been yet achieved in a 
wavefunction formalism, not even at large iV c . (See however the recent developments in 
Refs. [17,20].) On the other hand, the equations constructed in Refs. [4,9] propose such 
a synthesis at the level of the dipole scattering amplitudes (and for large N c ). 

Specifically, these equations refer to the scattering between a projectile which is a set of 
dipoles and a generic target. As shown in Ref. [11], the equations admit a natural Pomeron 
interpretation, where the "Pomeron" is the scattering amplitude for an external dipole: 
They describe the BFKL evolution of the Pomerons together with their interactions: 1 — > 2 
Pomeron splitting with the standard dipole vertex [22], and 2 — > 1 Pomeron merging with 
the "triple-Pomeron vertex" [33, 34]. By iterating these equations, one generates Pomeron 
loops [4, 11, 23]. So, from now on, we shall refer to them as 'the evolution equations with 
Pomeron loops'. 

The general structure of these equations alluded to above makes it tempting to search for 
a probabilistic interpretation in terms of a dipolar reaction-diffusion process, that is, a 
classical stochastic process in which a collection of dipoles — which represents either the 
target or the projectile — evolves through splitting and recombination. In what follows 
we shall briefly explain (leaving the details to the main body of the paper) why such an 
interpretation is tantalizing in the first place, what would be its benefits in practice, and 
why it does not really exist. But before we proceed, let us emphasize here that, although 
our main conclusion is essentially negative — we demonstrate the non-existence of a 
dipole interpretation — , this result has also positive implications, in the sense that it 
helps reorienting the efforts towards solving these equations (in particular, via numerical 
methods), and also towards better understanding the mechanism for gluon saturation. 
Besides, as we shall see, it helps clarifying some confusion which exists on this point in 
the literature. 

The reason why such an interpretation looks at the first sight appealing is because, taken 
separately, each of the terms appearing in these equations describes dipole splitting when 
considered in an appropriate Lorentz frame (and for large N c ). Specifically, the non- 
linear terms responsible for unitarity corrections can be seen as dipole splitting in the 
wavefunction of the projectile [5, 27, 28, 35], whereas the fluctuation terms relevant at low 
energy represent dipole splitting in the wavefunction of the target [3,4,22] (see Sect. 2 
below for details). Moreover, the Hamiltonian underlying the Pomeron loop equations, as 
constructed in Ref. [11], has a remarkable self-duality property [11,15], which in a given 
frame interchanges the splitting and the merging pieces of the Hamiltonian. Since the 
splitting piece has a well-founded dipolar interpretation, it is then tempting to search for 
a similar interpretation for the merging piece as well. 

Indeed, a fully dipolar picture of the evolution of the (projectile or target) wavefunc- 
tion, involving dipole splitting and merging, would be highly desirable, since very useful 
in practice. For instance, this would allow one to extend the Monte-Carlo method by 
Salam [26] up to arbitrarily high energy, and thus facilitate the numerical study of the 
evolution equations with Pomeron loops. More generally, this would establish an explicit 
correspondence between high-energy evolution in QCD and some Markovian "reaction- 
diffusion processes" which are intensively studied in the context of statistical physics (see 
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the textbook [36], and the review papers [37] for more recent developments), from where 
one could borrow some methods to QCD. Such a correspondence has been already proven 
to be useful, both at the level of the mean field approximation (the BK equation) [38], 
where it shed a new light on the 'geometric scaling' property of the solution [39-42], 
and in analyzing the role of the particle number fluctuations under asymptotic conditions 
(very high energy and arbitrarily small coupling constant) [3,4]. But, clearly, it would be 
extremely useful to extend this correspondence to a more interesting, non-asymptotic, 
regime, and rely on it in order to translate to QCD the huge experience accumulated with 
this problem in the context of statistical physics. Last but not least, a wavefunction pic- 
ture would allow for more general studies, and for the calculation of quantities other than 
the dipole scattering amplitudes (so like the gluon production in onium-onium scattering, 
a problem recently addressed in Ref. [43]). 

In view of the above, it may seem a little disappointing that such a dipolar picture 
cannot be actually given, for the reasons that we explain now. The problem comes from 
the gluon saturation effects, which cannot be described as dipole recombination. This 
problem can be seen already when trying to introduce saturation effects in the original 
dipole picture by Mueller [22,23]: when two dipoles within the same wavefunction are 
allowed to interact with each other by exchanging gluons, they produce not only dipoles, 
but also more complicated color states, like quadrupoles. This illustrates the fact that the 
merging processes responsible for saturation involve a non-trivial flow of color, and thus 
transcend the dipole picture. The actual theoretical tool for describing saturation, the 
JIMWLK equation [6-8], is written in terms of colorful degrees of freedom (gluons and 
the associated gauge fields), and cannot be directly simplified by taking the large- N c limit 
[11]. Such a simplification becomes possible only after applying the JIMWLK equation to 
dipole scattering amplitudes (thus yielding the Balitsky equations [5]); but in that case, 
the ensuing dipole interpretation refers to the evolution of the projectile through dipole 
splitting, and not to that of the target towards gluon saturation. 

This being said, the suggestive structure of the evolution equations with Pomeron loops 
alluded to above makes it tempting to search for an effective dipole picture at least, where 
by "effective" we mean a fictitious dipolar reaction-diffusion process which reproduces the 
actual equations for the scattering amplitudes without correctly reflecting the microscopic 
processes at work. (Such an effective picture would be, of course, still suitable for Monte- 
Carlo simulations.) And as a matter of facts, an effective dipole picture of this type 
has been recently proposed in Ref. [12], from the point of view of projectile evolution. 
However, as we shall demonstrate in what follows, that picture is truly ill-defined (it 
does not describe a physically acceptable reaction-diffusion process), and thus is useless 
in practice 4 . 

The basic problem is that the 2 — > 1 vertex which is supposed to describe dipole recom- 
bination within this effective picture (this is the same as the "triple-Pomeron vertex" 
[33,34]) has no definite sign, and thus it cannot be given a probabilistic interpretation, 



4 Besides, the analysis in Ref. [12] involves some errors of sign that we shall correct in our 
corresponding analysis in Sect. 5. 
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as a merging rate per unit rapidity. Very likely, this difficulty reflects the fact that the 
actual physical mechanism responsible for gluon saturation in QCD is a coherent, collec- 
tive, effect 5 — the suppression of the rate for gluon emission in the high density regime 
by strong color fields [8,31,44,45] — , and not a 'mechanical' reduction in the number 
of gluons. And indeed, the effective 'dipole recombination' process mentioned above has 
the property to leave the total number of dipoles in the system unchanged (unlike a real 
merging mechanism, which always decreases the number of particles). 

The above considerations motivate our strategy in this paper, which will be to demonstrate 
the non-existence of the dipolar picture through explicit construction. Namely, we shall 
consider a generic dipole-like reaction-diffusion process with arbitrary, non-local, vertices 
for dipole splitting and merging, and show that the only way for this process to reproduce 
the actual equations with Pomeron loops in QCD is by choosing the 'dipole recombination 
vertex' as the triple-Pomeron vertex mentioned before. Since the latter has no well-defined 
probabilistic interpretation, this choice renders the whole process meaningless. 

For our analysis to be convincing, we shall repeat it from the points of view of both target 
and projectile evolution. Indeed, a priori, these two points of view are not fully equivalent 
with each other, as they involve different 'factorization schemes' connecting the scattering 
amplitudes to the dipole densities: 

In the context of target evolution, it will be natural to evaluate the scattering amplitude 
for an external dipole in the two-gluon exchange approximation. This is motivated by the 
fact that (i) the dipole picture applies a priori to the dilute target regime, where this 
approximation is correct indeed, and (ii) when extended to the high-density regime, this 
approximation leads to a self-dual Hamiltonian [11], which has a suggestive, dipole-like, 
structure and thus is a natural starting point for the search of a dipole picture. 

From the perspective of projectile evolution, on the other hand, there is a priorino need for 
further approximations (except for the assumed reaction-diffusion process): Following a 
strategy by Kovchegov [27] (see also Refs. [12, 28,46]), the equations obeyed by the dipole 
densities can be exactly mapped onto equations for the scattering amplitudes, including 
multiple scattering effects to all orders. Yet, the two-gluon exchange approximation strikes 
back when one attempts to deduce the vertex for dipole recombination by comparison with 
the actual equations in QCD. 

Let us finally note that the lack of a probabilistic interpretation in terms of dipoles does 
not a priori preclude the existence of alternative stochastic descriptions, involving some 
more elementary degrees of freedom, so like the gauge field a a to be introduced in Sect. 
3. For instance, we have already mentioned that, within the JIMWLK formalism, the 
saturation processes have a well defined probabilistic interpretation [31], but in terms of 
gluon fields. Moreover, it is also known that the evolution equations with Pomeron loops 
can be compactly summarized into a single Langevin equation, with a rather unusual 



5 This is most transparently seen within the random walk formulation [8,31] of the JIMWLK 
evolution, where the probability for emitting new gluons at each step in the evolution involves 
Wilson lines and saturates to a field-independent value in the high-energy limit. 



5 



noise term though [9,10]; this is again suggestive of a classical stochastic process. An 
approximate version of this Langevin equation [4] has been recently used for numerical 
studies [47,48]. It would be very interesting to clarify the precise nature of the stochastic 
process underlying the equations with Pomeron loops, if any. 

The structure of this paper is as follows: In Sect. 2 we briefly review the evolution equations 
with Pomeron loops together with their physical interpretation from the perspectives of 
both target and projectile evolution. In Sect. 3, we study the possibility to interpret 
these equations in terms of dipole evolution in the target. To that aim, we rely on the 
Hamiltonian formulation introduced in Ref. [11] (to be reviewed in Sect. 3.1), that we 
first rewrite in terms of dipole operators (in Sect. 3.2) and then use to deduce evolution 
equations for the dipole densities (in Sect. 3.3). In the next section, these equations are 
compared to those generated by a genuine reaction-diffusion process. Specifically, in Sect. 
4 we introduce the Master equation defining a generic stochastic process with non-local 
vertices for splitting and recombination, and use it to deduce evolution equations for the 
dipole densities. (For convenience, the general equations are summarized in the Appendix.) 
Finally, in Sect. 5, we consider a reaction-diffusion process taking place at the level of the 
projectile. From the equations for the dipole densities (cf. Sect. 4), we deduce equations 
for the scattering amplitudes, which are eventually compared to the QCD equations with 
Pomeron loops. 



2 Evolution equations with Pomeron loops 

In this section, we shall briefly recall the evolution equations for dipole scattering ampli- 
tudes in QCD at high energy and large N c , together with their physical interpretation in 
terms of splitting and merging processes in either the target or the projectile [4,9, 10]. 

The relevant equations describe the rapidity evolution of the scattering amplitudes (T^) 
for the collision between a system of k dipoles (the 'projectile') and a generic hadronic 
target. The dipoles making up the projectile are assumed to be small enough for their 
scattering to lie within the scope of perturbation theory. In the course of the evolution, 
amplitudes (T^) corresponding to different values of k get mixed with each other, via the 
gluon splitting and recombination effects. Thus, the equations form an infinite hierarchy, 
whose general structure can be easily appreciated by inspection of the first two equations. 
These read [9] 



d(T(x,y)) 



27 



/ 



M xyz ® (T(x, y)) - M(x, y, z) (t*- 2 \x, z; z, y) 



) 



(2.1) 



dY 



z 



and, respectively, 
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d(T^ 2 \x 1 ,y 1 ;x 2 ,y 2 ] 
dY 



2tt 



M xiyiz ® (T (2) (a;i, y x ; x 2 , y 2 ) , 
-M(x 1 ,y 1 ,z) (T ( - 3 \x 1 ,z;z,y 1 ;x 2 ,y 2 )) + 

(Ct \ 2 O. I' 
2W 2^ J n ^y^ x ^y^ v ) ( 2 - 2 ) 



where a s = a s N c /ii, Y is the rapidity gap between the projectile and the target, {T(x, y)) 
is the scattering amplitude for an external quark-ant iquark dipole with the quark leg at 
x and the antiquark leg at y, {T^ 2 \x\, yi, x 2 , y 2 )) refers similarly to two external dipoles 
with coordinates (xi,yi) and (x 2 ,y 2 ), respectively, and so forth. 



The two kernels M and 1Z which enter the above equations are defined as follows: 

(x - y) 2 



M(x,y,z) 



(x - z) 2 (y - z) 2 



;2.3) 



is the dipole kernel [22, 24], which characterizes the differential probability for an elemen- 
tary dipole [x, y) to split into two dipoles (x, z) and (z, y) per unit rapidity. We have 
also introduced the shorthand notation 



M 



xyz 



f(x, y) = M(x, y, z)[-f(x, y) + f(x, z) + f(z, y)}. 



(2.4) 



Furthermore, TZ is the triple-Pomeron kernel (this name will be explained in Sect. 5.1) 
K.(x 1 ,y 1 ;x 2 ,y 2 \u,v) = J V 2 u V 2 v \M{u,v, z) A (x 1 ,y 1 \u, z) A (x 2 ,y 2 \z,v)~\, (2.5) 

z 

where a 2 Aq is (up to a sign) the amplitude for dipole-dipole scattering in the two-gluon 
exchange approximation and for large N c , with 



A (x,y\u,v) = I 



In 



(x — v) 2 (y — u) 21 2 
(x — u) 2 (y — v) 2 



(2.6) 



The general structure of the higher equations in this hierarchy can now be easily under- 
stood : The amplitude for k dipoles with k > 2 is coupled to the one for k— 1 dipoles, to 
that for k dipoles, and to that for dipoles. For more clarity, in what follows we shall 
refer to these three types of terms as the fluctuation terms, the (normal) BFKL terms, and 
the unitarity corrections, respectively. (These names should become clear from the subse- 
quent discussion.) If one ignores the fluctuation terms, which are formally suppressed by 
an additional factor of a 2 (see, e.g., the last term in Eq. (2.2)), then the above hierarchy 
boils down to the large-iV c version of the Balitsky-JIMWLK hierarchy [5,28,35]. How- 
ever, as pointed out in Refs. [2-4], the fluctuation terms are in fact leading-order effects 
whenever the target is so dilute (or, equivalently, the external dipole (x, y) is so small) 
that (T{x, y)) < a 2 . Note also that T = 1 is a fixed point of the evolution described by 
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these equations 6 , which physically corresponds to the 'black disk' limit at high energy. 



The above equations are invariant under a Lorentz boost along the collision axis [11], 
but the physical interpretation of the various terms depends upon the frame in which we 
visualize the evolution. As mentioned in the Introduction, it is always possible to choose 
the frame in such a way that a given term corresponds to dipole splitting. 

Namely, the BFKL terms and the unitarity corrections (that is, all the terms in Eqs. (2.1)- 
(2.2) except for the last term in Eq. (2.2)) can be interpreted as dipole splitting in the 
evolution of the projectile : If the rapidity increment dY is given to the projectile, one of 
the k original dipoles there splits into two new dipoles, which then scatter off the target. 
From this perspective, the BFKL terms ~ correspond to the separate scattering of 
the daughter dipoles (plus virtual corrections), whereas the unitarity effects ~ are 
rather associated with their simultaneous scattering. 

Alternatively, the rapidity increment dY can be given to the target, in which case it is 
the fluctuation term ~ T^ -1 ) (e.g., the last term in Eq. (2.2)) which admits a simple, 
dipolar, interpretation [4]: In the weak scattering regime where such fluctuations are 
truly important, the target is dilute and can be itself described, within the large-iV c 
approximation, as a collection of dipoles which evolve through dipole splitting [10, 22, 24]. 
Then the last term in Eq. (2.2) describes the splitting of the target dipole (u, v) into 
two new dipoles (u, z) and (z,v), which then scatter off the external dipoles («Ci,j/i) 
and (#2, 2/2) Yia - twice two-gluon exchange. This interpretation becomes perhaps more 
transparent after performing some integrations by parts to rewrite this term as 



(this is the original form of this term in Ref. [9]) and then observing that, in the dilute 
regime, V^Vj, (T(u, v)) is proportional to the average dipole number density in the target 
(n(u,v)), as it will be explained in Sect. 3.2 below. 

The BFKL terms too can be understood as dipole splitting in the target [22,24]. On the 
other hand, the unitarity corrections do not admit a dipolar interpretation in this frame: 
From the perspective of target evolution, they rather correspond to saturation effects, that 
is, to gluon recombination in the target wavefunction. Similarly, the fluctuation terms can 
be associated with gluon recombination in the projectile wavefunction. We thus arrive at a 
physical picture in which, in a given frame, some of the terms in the evolution equations at 
large N c can be directly interpreted as dipole splitting, whereas the other terms correspond 
(in a less direct way, though) to gluon merging. As anticipated in the Introduction, and it 
will be explained in detail in what follows, the gluon merging is the process which prevents 
a fully dipolar interpretation for the evolution equations even at large N c . 



6 This is perhaps more obvious after rewriting the fluctuation terms as in Eq. (2.7). 




M(u, v, z) A(#i, V\\u, z) A (x 2 , y 2 \z, v) 



x V 



y 2 v (t(u,v)), 



(2.7) 
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3 Dipoles in the target 



In this section, we shall demonstrate that Eqs. (2.1)-(2.2) for the scattering amplitudes 
cannot be obtained as the result of a stochastic process involving dipoles in the target. To 
that aim, we shall start with the Hamiltonian formulation of these equations [11], which 
in the dilute regime reduces to the original dipole picture by Mueller (in its 'color glass' 
version of Refs. [10,24]), and thus represents a natural starting point in the search for a 
dipole interpretation. 



3.1 Hamiltonian formulation in the a -representation 



As shown in [11], the whole hierarchy of evolution equations for dipole scattering am- 
plitudes as presented in Sect. 2 can be derived from the following effective Hamiltonian 



H = H + H^ 2 + # 2 ^i, (3.1) 

where the three distinct terms correspond, respectively, to BFKL evolution, gluon split- 
ting, and gluon merging in the wavefunction of the target. They are given by 7 : 

"» = ^?i / M <»■ '"•<"> - Q °w '"•('> - Q >" J ( u) J(v) • (3 ' 2) 



#1- 



16iV c 3 



^ j M{u,v, z)g(u 1 \u, z)Q(v 1 \u, z)g(u 2 \z,v)g(v 2 \z,v) 

72T72~c(„\„.c(..\ 6 6 6 6 

Sa a (ui) Sa a (vi) 5a b (u 2 ) 5a b (v2) ' 



xV 2 u V 2 v a c (u)a c (v) 



(3.3) 



c u,v,z v ' v ' 

In these equations, a a (z) is the total color gauge field radiated by color sources within the 
target (in an appropriate gauge), and the function Q(ui\u,z) represents (up to a factor 
gt a ) the classical field created at ui by the elementary dipole (u, z): 

G{ui\u,z) = j-ln I" 1 ~ (3.5) 
An (ui — u) 2 

In what follows, we shall refer to the above form of the effective theory as the "en- 
representation". The operators a a (z) and 5/5a a (z) can be interpreted as 'Fock-space 

7 The splitting piece Hi— ^ of this Hamiltonian, Eq. (3.3), has been first obtained by Mueller, 
Shoshi and Wong [10]. The merging piece H 2 ^i, Eq. (3.4), has been introduced in Ref. [11], as 
the dual counterpart of H\^2- 
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operators' which annihilate and create a gluon, respectively. This makes it clear that the 
above Hamiltonian describes the evolution of the gluon distribution in the target, which 
proceeds through BFKL evolution (H ), 2^4 gluon splitting (if 1 _ +2 ), and 4^2 gluon 
merging (if 2 ^i)- From this perspective, the splitting Hamiltonian (3.3) [10] should be 
more naturally denoted as i^2->4 (and similarly H^ 2 for the merging term (3.4)), but the 
present notations anticipate the 'dipolar' form of the Hamiltonian, to be presented later. 
The Hamiltonian in Eqs. (3.1)-(3.4) has a remarkable self-duality property [11, 15], to be 
discussed in Sect. 5.1. 



In this representation, the scattering amplitudes for external dipoles are obtained as 
(TW) = (T(1)T(2)...T(k)), where T = T(x,y) is the following scattering operator: 



which appears to describe single scattering in the two-gluon exchange approximation. 
However, within the present formalism, Eq. (3.6) is to be used in the whole kinematical 
range, including at high energy where the multiple scattering is important. This is ap- 
propriate since the a-representation described here is only an effective theory, in which 
the various microscopic process are not always faithfully depicted, but which reproduces, 
by construction, the correct evolution equations for the scattering amplitudes (that is, 
the hierarchy starting with Eqs. (2.1)-(2.2)). Indeed, by using the above formulae, it is 
straightforward to check that the Hamiltonian evolution equations 



generate the expected hierarchy, so long as one neglects contractions which give terms 
suppressed at large- N c . For example, when acting with S 2 / Sa^da^ from Hi^ 2 on the 
two-dipole scattering operator T^(xi,yi;x2,y 2 ), one should keep in the equation only 
the terms which are generated when both derivatives act on the same dipole. In that sense, 
the action of the Hamiltonian is constrained by a large-iV c counting rule. Moreover, the 
Hilbert space for the a-representation is also constrained: Because of the assumptions 
underlying the construction of the Hamiltonian (3.1)-(3.4) (see the discussion in Ref. 
[11]), this Hamiltonian cannot be used to act on arbitrary operators built with a, but only 
on the operators built with the scattering amplitude (3.6). (Such acceptable operators are 
gauge invariant in the sense explained in Ref. [49].) 

3.2 Hamiltonian theory in the D -representation 

In the a-representation introduced above, the dipolar structure of the Hamiltonian is not 
explicit. But as we shall explain below, this structure becomes manifest after a change 
of representation, which has also the advantage to eliminate all the constraints on the 
Hamiltonian structure that we have previously mentioned. 

Specifically, one can reexpress the Hamiltonian and the interesting observables (the scat- 




(3.6) 




(3.7) 
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tering amplitudes) in terms of operators which are bilinear in the number of fields a or of 
functional derivatives 8 /8a. If p a denotes the color charge density in the target, so that 



-V 2 x a a {x) = p a (x), 

then we define the D-representation via the following identifications: 
1 



D x y 



g 2 N c 



pa(x)p a (y) 



x^y 



(3.8) 



(3.9) 



D 



xy 



1 + 



4iV r 



8 



8p a (x) 8p a (y) 



(3.10) 



The notation is slightly abusive, since D and are not Hermitian conjugate to each 
other in any sense, but it is suggestive of the fact that these operators act as creation and 
annihilation operators for dipoles in the target (see Sect. 3.3 below). 

When expressed in terms of D and D\ the Hamiltonian takes a rather compact form: 



H 



2^ 



2 V2tt/ 2tt 



M(u,v,z)D uv [-Dl v + Dl z Dl 

/ n(u 1 ,v 1 ;u 2 ,v 2 \u,v) D UlVl D 

U2V2 



(3.11) 



where the expression in the first line has been obtained by summing up the BFKL 8 and 
the splitting pieces of the Hamiltonian, cf. Eqs. (3.2)-(3.3), while that in the second line 
is a rewriting of the merging Hamiltonian, Eq. (3.4). 

Furthermore, the dipole scattering amplitude can be related to the operator D xy via 
T(x,y)[D] = a] J A (x,y\u,v)D uv , (3.12) 

which follows from Eqs. (3.6), (3.8), and (2.6), together with the condition that the system 
be globally color neutral: J u p a (u) = 0. 

To complete the definition of the new representation, one still needs to specify the algebra 
of the operators D and DK To that aim, we rely on their respective expressions in terms 
of p and 8/8p, cf. Eqs. (3.9)-(3.10), to conclude that, in the representation in which D xy is 
diagonal (i.e., a pure function) and in the large- N c approximation, can be represented 
as a functional derivative with respect to D : 



with 



xy 



8D 



xy 



D uv — ^ (8xu8yv ~\~ 8 X v8yzi) ■ 



(3.13) 



There is in fact a slight subtlety concerning the BFKL Hamiltonian, which in going from 
Eq. (3.2) to Eq. (3.11) has been rewritten in a different, but equivalent, form; see the discussion 
in Ref. [11]. 
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To understand how this arises, note first that, for u ^ v, 



9 



2 



4N r 



1 2 



5p a {x) 5p a (y) 



1 



Pa(u)p a (v) ~ - (5 ), (3.14) 



g 2 N c ' ay /,ay ' " 2 



where we have used N% — 1 ~ N% at large N c . Furthermore, if we temporarily return 
to the p-representation in Eqs. (3.9)-(3.10) and consider the action of on a string of 
.D-operators, so like D\. y D ulvl D U2V2 . . . , then it is clear that the dominant result at large 
N c is the same as obtained when acting with the functional derivative 5/5D xy . 

Thus, by using the expression (3.11) for the Hamiltonian together with Eq. (3.13) for D\ 
one ensures that the action of H on operators built with D generates only the would-be 
dominant contributions at large- N c . Moreover, the Hilbert space in the D-representation 
is not constrained anymore: Indeed, Eq. (3.12) can be inverted to yield 

D*y = - 4 V 2 x V 2 y T(x,y) for x + y, (3.15) 

showing that any operator built with D is physically acceptable. Thus, as anticipated, by 
changing from the a to the D-representation we have eliminated all the constraints. 

It is now a straightforward exercise to verify that the evolution equations for (D^) = 
(£>(1)£>(2) . . . D(k)) generated as: 

d{D(K)) = (H £><">) (3.16) 



dY 

lead to the expected equations for the scattering amplitudes after also using Eq. (3.12). 



3.3 A dipolar picture and its limitations 



We now turn to the dipolar interpretation of the evolution described by the effective 
Hamiltonian (3.11). As we shall soon discover, such an interpretation is truly meaningful 
(in the sense of having a probabilistic meaning) only in the dilute regime, where it reduces 
to the standard dipole picture by Mueller [22,23]. Namely, for a dilute target, the recom- 
bination effects can be ignored and the Hamiltonian (3.11) reduces to the expression in its 
first line, i.e., i^dipoie = H + H^ 2 , which is recognized as the 'color glass' version of the 
Hamiltonian generating the dipole picture [10,24]. As we briefly recall now, the evolution 
driven by i^dipoie leads to a probabilistic picture of the target, as a collection of dipoles 
which evolves through splitting (see Refs. [10, 19,22-24] for more details). 

Specifically, in the dilute regime, the operators and D introduced before can be in- 
terpreted as dipole creation and annihilation operators, respectively, with the following 
commutation relation (at large N c ): 



D 



($xu$yv + $xv$yu) ■ (3-17) 
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The target distribution is then obtained by successively acting with the creation opera- 
tors on the 'dipole vacuum', and can be represented as the following 'color glass weight 
function' : 

oo . N 

MP] = E dr N P N ({u i ,v i };Y)l[D\u i ,v i )5[p}, (3.18) 

JV=1 i=l 

where the functional distribution 5[p\ represents the vacuum (no color charge, and hence no 
dipole !), Z]y({ui, Vi}) = ]\f =1 D^ [ui,Vi)5[p] describes a given configuration of N dipoles, 
Pn({u{, Vi}; Y) is the corresponding probability at rapidity Y, and dT^ denotes the mea- 
sure for the iV-dipole phase-space: dTjy = (1/N\) Ylf =1 d 2 Uid 2 Vi. 

The 'weight function' (3.18) acts as a probability distribution for p a when computing 
target expectation values. For instance: 

(D(x,y)) = J D[p]D(x,y)Z Y [p] 

oo N 

= E / dT N P N {{ Ui , Vi};Y)Y,KJv iy = (n(x,y)), (3.19) 

N=l J i=l 

where in writing the second line we have used Eqs. (3.17) and (3.18), and then recognized 
the definition of the average dipole number density. 

In view of the discussion above, the dipole interpretation of the evolution generated by 
-f^dipoic is quite transparent (cf. the first line of Eq. (3.11)) : The term involving D^D^D^ 
describes the splitting of the original dipole (u, v) into two new dipoles (u, z) and (z, v), 
with a probability density (a s /2ir)JA(u, v, z) per unit of phase-space and per unit rapid- 
ity. The other term, proportional to —D UV D^ UV , describes the decrease in the probability 
that the original dipole remain unchanged after one step in the evolution. 

To verify that this probabilistic interpretation is mathematically well founded, one can use 
-f^dipoie to deduce evolution equations for either the probability densities Pjv, or the dipole 
number densities (n^) = (D^), and then check that the ensuing equations have the 
structure expected for the Master equations associated with a stochastic process. For the 
splitting problem described by iJdipoie, both tests have been performed in the literature 
[4,24,26], and they confirm the probabilistic interpretation given above. In what follows 
we shall follow the second strategy, namely, the one which focuses on the particle densities, 
but we shall do so in the more general context of an evolution involving both splitting and 
merging, as generated by the full Hamiltonian in Eq. (3.11). 

Indeed, a question which naturally arises at this point is whether one can extend this dipole 
picture to the high-density regime as well, that is, to the merging piece of the Hamilto- 
nian, as given in the second line of Eq. (3.11). Since proportional to D uxvx D u ^ m D^ uv , 
is seems natural to interpret this piece as describing the recombination of two dipoles 
(ui,Vi) and (1*2,^2) into a single dipole (u,v), with a probability density which is pro- 
portional to 1Z(ui, V\\ U2, V2\u, v) (and with an overall sign which is not clear yet; this will 
be clarified in Sect. 4.2). However, a more careful analysis reveals that this interpretation 
is in fact deceiving, as shown by the following argument: For the probabilistic interpre- 
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tation to be meaningful, the 'triple-Pomeron' kernel TZ must have a definite sign (since a 
probability density cannot change sign !) However, a brief inspection of Eq. (2.5) reveals 
that the actual kernel 1Z has not a fixed sign. Indeed, using the fact that this kernel is a 
total derivative and that there are no boundary terms when it is integrated over all the 
transverse space (as we shall shortly check), we immediately find 

J 7£(wi, vi, u 2 ,v 2 \u, v) = 0, (3.20) 

showing that the value of the vertex must be negative in some regions of the transverse 
space, and positive in some others. Clearly, this property prohibits any probabilistic in- 
terpretation for the merging vertex in the Hamiltonian (3.11). 

At this point, we should stress that the aforementioned sign property of the triple- 
Pomeron kernel does not entail any sign problem for the fluctuation term in Eq. (2.2), and 
thus for the actual evolution equations in QCD: Even though 1Z has no definite sign, the 
convolution yielding the last term in Eq. (2.2) is nevertheless positive-definite within the 
whole kinematical range in which this term is important, that is, in the weak scattering 
regime. To verify this, it is preferable to rewrite this term as in Eq. (2.7) and then use 
the fact, for a dilute target, the quantity V^V^ (T(u,v)) is positive-semidefinite, since 
proportional to the average dipole number density (cf. Eqs. (3.15) and (3.19)). 

To render the previous arguments more precise, we shall compare in Sect. 4 the equations 
for the dipole number densities generated by the Hamiltonian (3.11) to those arising from 
an actual reaction-diffusion process with generic vertices for splitting and recombination. 
To that aim, let us derive here the corresponding equations in QCD. The 'dipole number 
density' is defined as the expectation value of D, cf. Eq. (3.19), and similarly for the 
higher /t-body densities (e.g., (DD) is the average pair density, etc.). One should keep in 
mind, though, that such expectation values can be meaningfully interpreted as particle 
number densities only in the dilute regime. The corresponding equations are then obtained 
according to Eq. (3.16). It is convenient to separate the respective contributions of the 
splitting (H + Hi^ 2 ) an d merging (if 2 ->i) terms in the Hamiltonian. 

For the BFKL+splitting contributions, one finds: 



d{D x 



dY 



l->2 



2^ 



IC UVZ ® (D uv ) , (3.21) 



and furthermore 



dY 



l->2 27T 



J fcuiviz ® (D UlVl D 



U2V2 I 



+5(u 2 - Vi)M(ui, v 2 , u 2 ) (D. 



U1V2 I 



+ {l <-> 2}, (3.22) 



where he have used the shorthand notation 



1C*vz ® /(«, v ) = -M(u, v, z)f(u, v) +M(u, z, v)f(u, z) 

+M(z,v,u)f(z,v). (3.23) 
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Note, in particular, the term linear in (D) in the r.h.s. of Eq. (3.22) for (D^). This term 
describes a 'fluctuation' in which two dipoles which have a common leg are produced in 
one step of the evolution, via the splitting of a single original dipole. Through Eq. (3.12), 
this mechanism produces the fluctuation term in Eq. (2.2). 



The recombination piece in the equation for (D) is similarly obtained as 

J TZ(u 1 ,v 1 ;u 2 ,v 2 \u,v) (D U1V1 D U2V2 ) . (3.24) 



d(D uv ) _ ( a£\ 2 cts_ 
2tt) 2ir 



dY 



2->l 



Via Eq. (3.12), this piece produces the non-linear term ~ (T^) in the r.h.s. of Eq. (2.1). 
The higher equations in the hierarchy for (D^) can be derived similarly, but the equations 
written above will in fact suffice to illustrate our point in the subsequent discussion. 

To conclude this section, let us explicitly check the fact that the surface terms for the inte- 
gral (3.20) vanish indeed, as announced. Concentrating, for example, in the u integration 
there, we have 

j lZ( Ul , «i; « 2 , v 2 \u, v) = j V 2 u F(u) = Jiim j dS- V u F(u), (3.25) 

U u C R 

with the one-dimensional "surface" integral to be taken along Cr, the circle of radius R. 
More precisely 

dE = wd0u, (3.26) 

where is the polar angle of the vector u, u is the unit vector along the radial direction, 
and u = |u| is equal to R on Cr. By expanding the dipole kernel M.(u,v,z) and the 
dipole-dipole scattering Aq{ui, V\\u, z) for large values of u it is straightforward to show 
that 

F(u) = A + + . . . =>- V u F(u) = (cos0u + sin00) + ..., (3.27) 



u u 



with the scalar A and the vector B being independent of u. Since we are going to integrate 
over 0, we have assumed, without any loss of generality, that the polar angle of the vector 
B is zero. Now it is trivial to see that the insertion of Eqs. (3.26) and (3.27) into Eq. (3.25) 
gives a result which vanishes indeed when R — > oo : 

27T 

/ TZ(u\, V\, u 2 , v 2 \u, v) — — B I d(f)COS(f) lim — = 0. (3.28) 

J J R^oo R 



4 The general reaction diffusion process 



In this section, we shall describe a generic reaction-diffusion process with non-local ver- 
tices for particle ('dipole') splitting and recombination. Starting from the Master equation 
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for probabilities, we shall deduce an hierarchy of evolution equations for the particle num- 
ber density and its correlations. For more clarity, we shall discuss separately the effects 
of splitting and merging, which is indeed possible since the two types of processes take 
place independently from each other. (In the complete equations, which are displayed for 
convenience in the Appendix, the corresponding terms simply add up to each other.) This 
is in line with our general strategy which is to emphasize the disparity (in so far as the 
dipolar interpretation is concerned) between splitting and merging contributions to the 
actual evolution equations in QCD. 



4-1 Master equation for splitting 



To facilitate the correspondence with the QCD problem at hand, we shall use a QCD- 
inspired terminology also for the generic reaction-diffusion problem. The system is there- 
fore a 'hadron' whose general state is characterized by the number of 'dipoles' N together 
with the associated coordinates {u\ ... un), where U{ is a 4-dimensional variable represent- 
ing in a collective way the coordinates of the quark and the antiquark of the i-th dipole. 
Now, let us assume that, under a step dY in 'time' (or 'rapidity'), a dipole can split into 
two non-locally. The corresponding Master equation is 



dP, 



N 



dY 



dP, 



(<?) 



N 



dY 



+ 



dP%> 



dY 



(4.1) 



where P/v = Pn(ui---Un) is the probability density for a given configuration of iV dipoles 
(this is totally symmetric under the permutations of the dipoles) and the superscripts (g) 
and (£) stand for the gain and loss terms, respectively. In their most generic form, the 
two terms on the r.h.s. of Eq. (4.1) are given by 



dP { N 9 \ Ul . 



■u N ) 



dY 



Y J a(v\uiUj) P N - 1 (ui...^ i ...^ j ...u N v), 

i^j v 



(4.2) 



and, respectively, 
dP$(u 1 ...u N ) 



dY 



Y J a(ui\vw) P N (ui...u N ), 



(4.3) 



with a(v\uiUj) = a(v\ujUi) the differential probability that the dipole v split into the 
dipoles Ui and Uj (this is necessarily positive semi-definite, and symmetric in the coor- 
dinates of the child dipoles), and where it is understood that a slashed variable is to be 
excluded. It is straightforward to show that Eqs. (4.1)-(4.3) conserve the total probability, 
i.e. 



_d_ 

dY 



00 . 

Y / dT N P N (ui...u N ) = 0, 

N=l 



(4.4) 



where the differential space of integration is 
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&T N = —du l ..Au N = —. 



(4.5) 



(The factor 1/N\ is needed to divide out the number of equivalent configurations). In 
fact the probability conservation property is more stringent, as it holds already before 
summing over N : 



/ 



dr 



dp\ e) 



N 



N 



dY 



+ 



ap(9) 

dr- dF " +1 



N+l 



dY 



0. 



(4.6) 



Using the above Master equation, we would like to derive the evolution equation for a 
generic operator O which depends upon the transverse coordinates of the dipoles. The 
expectation value of such an operator may be written as 

oo „ 

(O) = E / dr N P N ( Ul ...u N ) O n ( Ui ...u n ). (4.7) 

N=l J 

By using the Master equation together with the assumed permutation symmetries for Pn 
and for the splitting vertex, we find 



d(0) 



dY 



E /^--ni 0((u N \u N+1 u N+2 )P N (u 1 ...u N ) 

N=l 1 J- 



X 



O n (u 1 ...u n ) + O n+1 (u 1 ...u N - 1 u n+1 u n+2 ) ■ (4.8) 



To be more specific let us first write the equation which describes the evolution of the 
average dipole density. According to the general definition (4.7), this is defined as 

oo „ 

(n(v)) = E dr N P N (u 1 ...u N )n N (u 1 ...u N ), (4.9) 

N=l J 

where un is specified to be 

N 

n N = Y / S(u l -v). (4.10) 
i=\ 

Then it is easy to see that only three of the (^-functions survive in the square bracket in 
Eq. (4.8). These give rise to three terms in the equation for (n(v)) and we arrive at 



d {n(v)) 



dY 



J ^2a(w 1 \vw 2 ) - a(v\w 1 w 2 ) (n(v)) , (4.11) 



WlW'2 



where the factor of 2 in the first term comes from the symmetry property of the vertex 
a{v\wiW 2 ) = a(v\w 2 Wi). The interpretation of the r.h.s. of Eq. (4.11) is transparent: there 
is a positive contribution, since the arbitrary dipole W\ can split into the dipoles v and 
w 2 , but the latter is not "measured" though, and there is a negative contribution, since 
the dipole v can be lost by splitting into w± and w 2 . 

Similarly, the dipole-pair density is defined through 
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n 



(2) 
N 



= J25(u t -v 1 )5(u 3 -v 2 ), (4.12) 



and evolves according to 



dY 



= / \2a(w 1 \v 1 w 2 ) (n ( - 2 \w 1 v 2 )\ - a(v 1 \w 1 w 2 ) (V 2) (^2) 

+ {vi ^v 2 ] + 2 J a(w\v lV2 ) (n(w)) . (4.13) 



All the terms in the r.h.s. except for the last one are the straightforward generalization 
of the corresponding terms in Eq. (4.11). As for the last term in Eq. (4.13), this describes 
a 'fluctuation' in which an arbitrary dipole w splits into two dipoles v\ and v 2 which are 
both measured. 

It is now easy to understand the general structure of the hierarchy corresponding to 
splitting: The r.h.s. of the equation for (n^) (the K-dipole density, defined by analogy 
to Eq. (4.12)) involves terms proportional to (n( K )), associated to splitting processes in 
which only one of the child dipoles is measured, and terms proportional to (n*^ 1 )), in 
which both child dipoles are measured. The latter are the 'fluctuation terms'. 

It is furthermore straightforward to check that the above equations for (n) and (n^) will 
reduce to Eqs. (3.21) and (3.22) of the QCD problem for the specific vertex 

a{v\wiw 2 ) = «3p(f v\wiWi, w 2 w 2 ) 

= I 7T 6 ( v ~ w i) s (v - W2)S(w! - w 2 )M{v, v, tox) + f 1 <-> 2} , (4.14) 

(the subscript "3P" will be explained in Sect. 5.1) where we have separated each 4- 
dimensional variable into two 2-dimensional ones, denoting the coordinates of the quark 
and the antiquark, respectively. Therefore, the splitting-in-the-target processes in QCD, 
as described by the first line in the effective Hamiltonian (3.11), have an equivalent de- 
scription in terms of a Master equation, and thus a well-defined (dipolar) probabilistic 
interpretation, as anticipated. The same holds therefore for the BFKL+fluctuation terms 
in the equations for the scattering amplitudes, cf. Eqs. (2.1)-(2.2), when these are viewed 
from the perspective of target evolution. 



4-2 Master equation for recombination 



Let us now similarly consider the master equation for a 'dipole' system which evolves only 
through recombination. In analogy with Eqs. (4.1)-(4.3) we can write 



dP, 



N 



dY 



dP 



(9) 



N 



dY 







+ 



dY 



(4.15) 



where the gain and loss terms have now the following structure 
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dP { N 9 \ Ul ...u N ) 



dY 



E 

i , 



/3(vw\u i )P N+1 (u 1 ...^ i ...u N vw), 



(4.16) 



dP$(u 1 ...u N ) 



dY 



= - E/ P{u i u j \v)P N {u 1 ...u N ), 



(4.17) 



where f3{uiUj\v) = f3(ujUi\v) is the (positive-semidefinite) differential probability that 
the dipoles Ui and Uj recombine with each other and thus produce the dipole v. The 
probability conservation law (4.4) is again satisfied, and the analog of Eq. (4.6) reads 



/ 



dr 



dPl 9) 



N~ 



N 



dY 



+ 



t dP (e) 

jdr - (,lx - 



N+r 



dY 



0. 



(4.18) 



The evolution equation for the generic operator O is obtained as 



d(0) 



dY 



E 



dr 



7V+2 



P (U N U N+1 1 U N+2 )P N («! . . .U N+ i ) 



X 



(4.19) 



O n+ i(ui...u n+ i) + O n (ui...u N -iu n+2 ) 
After straightforward algebra we arrive at the following equation for the average density 
d(n(v)) _ f r na , l.. \ /„(2) i N \ , at „.. i „.\ /„ (2) i Al (42Q) 



dY 



= J ~2P(wiv\w2) (n^ 2 \wiv)^ + P(wiW 2 \v) (n (,2 \wiW 2 ) 



and for the average pair density 

d (n <y2 \viV2) 



dY 



\—2f3(w-iV\\w2) (n^ (wiviv 2 )} + P(wiw 2 \vi) (n^ '(wiw 2 v 2 ) 
+ {vi <-> v 2 ) - 2 J (3{ Vl v 2 \w) (n^{ Vl v 2 )) . (4.21 



As manifest on these equations, the recombination process has two opposite consequences: 
on one hand, it yields a negative contribution to, say, the average number density (n(v)), 
because a dipole v can merge with any other dipole W\ which is available in the system, 
and thus disappear; on the other hand, it yields a positive contribution to (n(v)), because 
two arbitrary dipoles w 1 and w 2 can fuse with each other and thus create the dipole 
v. Besides, for higher-point correlations (n^) with k >2 there are additional negative 
contributions since any two among the k external dipoles can recombine with each other. 

Yet, it is easy to check that the negative contribution to d{n(v))/dY prevails after in- 
tegration over v , and thus the global effect of the recombination is to reduce the total 
number of dipoles in the system, as intuitive. Indeed, Eq. (4.20) implies: 



/ 



d {n(v)) 



dY 



(3(wiw 2 \v) (rS 2 ^ 



\WiW 2i 



< 0. 



(4.22) 



vwiW2 
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(By a similar argument, one can convince oneself that the global effect of splitting is to 
increase the total particle number.) In particular, if the vertex is fully local, /3(wiW2\v) oc 
5(wi —v)5(w2 —v), then the effect of the recombination is to decrease the particle number 
density at any v. 

Now let us compare the above equation for (n(v)) to the QCD equation (3.24), as gen- 
erated by the merging piece (the second line) in the Hamiltonian (3.11). One then finds 
that, formally, Eq. (3.24) can be written in the form (4.20) by choosing (3(w 1 w 2 \v) = 
/3 3 f>(w 1 w 1 ;w 2 w 2 \vv) with 

P 3W >(u 1 Vi;u 2 v 2 \uv) = (J^J ^Tl{ui,vi; u 2 ,v 2 \u, v) + {l ^ 2} . (4.23) 

Indeed, with this choice, the first term in Eq. (4.20) vanishes identically (recall that 1Z is 
a total derivative w.r.t. its arguments u and v), while the second term there becomes the 
same as the r.h.s. of Eq. (3.24). However, it is obvious that such a choice for (3 is physically 
meaningless: (i) The vertex (4.23) has no fixed sign, as already discussed after Eq. (3.24), 
and thus cannot be interpreted as a merging rate, (ii) With this choice, the would-be 
dominant term in the r.h.s. of Eq. (4.20) (namely, the term with an overall factor —2) 
appears to vanish, and the whole non-trivial contribution comes from the second term 
there, (iii) For f3 = flsp, the r.h.s. of Eq. (4.22) is exactly zero, showing that the QCD 
equation (3.24) cannot be really interpreted as describing dipole merging. 

We thus conclude that, when the evolution is viewed as gluon evolution in the target, the 
non-linear terms (unitarity corrections) in the equations for the scattering amplitudes 
cannot be attributed to a process like the recombination of dipoles. In the following 
section we shall see that a similar problem appears also when we evolve the projectile, 
except that, in that case, the problem arises for the fluctuation terms. 



5 Dipoles in the projectile 

In this section we shall consider Eqs. (2.1)-(2.2) for the scattering amplitudes from the 
point of view of projectile evolution and show that, in this context too, the dipole inter- 
pretation is plagued with the difficulty identified before: namely, it necessarily involves 
the 'dipole recombination' vertex of Eq. (4.23), which has no probabilistic meaning. We 
shall discuss two a priori different strategies to introduce dipoles in the projectile, and 
show that they both suffer from the deficiency alluded to above: 

First, in Sect. 5.1, we shall reformulate the Hamiltonian underlying the QCD equations (cf. 
Sect. 3.1) in terms of new degrees of freedom which are naturally interpreted as dipoles 
in the projectile. Then, in Sect. 5.2 we shall assume — following Ref. [12] - that the 
projectile evolves according to a dipolar reaction-diffusion process, which involves both 
splitting and recombination. The equations for the scattering amplitudes that we shall 
obtain in this way will be finally compared to the actual QCD equations in Sect. 5.3. 



20 



5.1 Hamiltonian in the S -representation 



Returning to the evolution Hamiltonian introduced in Sect. 3.1, we shall rewrite it here 
in a different form — the "S'-representation" — , in which the explicit degrees of freedom 
are dipoles in the projectile, or, more precisely, the ^-matrices for the scattering between 
such dipoles and the target. Specifically, we introduce the following operators 



,2 



SL = 1 " iw [«"(*) " « a (2/)f (5-1) 



which in the large-iV c limit satisfy 

Sxyi 'S'ttii ~ ($xu$yv H~ $xv$yu) i (5-3) 

and in terms of which the Hamiltonian reads 
H 



u,v,z 



s ^] = ^ J M(u,v,z)[-Sl v + Sl z Sl 

J K{u u t>i; u 2 , v 2 \u, v) S j uv S UlVl S U2V2 , (5.4) 



1 / a 8 \ 2 a s 
~2 Vhr) 2^ 

where the expression in the first line is the sum Hq + i?2-»i, cf. Eqs. (3.2) and (3.4), while 
that in the second line is a rewriting of Hi^ 2 , Eq. (3.3). 

Note that the above S-form of the Hamiltonian is dualto its previous D-form, Eq. (3.11), 
where the duality transformation [11,15] consists in replacing 

D^S , — > S\ (5.5) 

and then exchanging the relative ordering of the operators 9 5' and S. Note also that 
the splitting and merging pieces in the original Hamiltonian (3.2)-(3.4) get interchanged 
by this transformation: The merging piece H 2 ^i, which is the same as the second line of 
Eq. (3.11), is dual to the expression in the second line of Eq. (5.4), which represents the 
splitting piece Hi^ 2 . This symmetry reflects the s elf- duality of the complete Hamiltonian 
under the following transformations [11]: 

5 1 5 

^ ~ 'Jpjxj ' 7 I^x) " Pa{x) ' (5 - 6) 

(followed by Hermitian conjugation), a property that has been actually used in Ref. [11] 
to construct the effective Hamiltonian. As it should become clear from the subsequent 



9 This permutation corresponds to Hermitian conjugation, in the sense of an integration by parts 
in the color glass functional integral yielding the expectation value over the target wavefunction 
(see the discussion in Ref. [11]). 
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discussion, the physical meaning of the duality transformation is to exchange dipoles in 
the target with those in the projectile. 

Indeed, the operator is recognized as the S'-matrix for the scattering between the 
projectile dipole (x, y) and the target in the two-gluon exchange approximation. Since 
the average amplitude (T) is obtained as (T) = (1 — S^), it is in fact preferable to work 
in the S^-representation, in which is a pure function to be denoted as S, and where 
T = 1—S. This is what we shall properly mean by the "S'-representation" in what follows. 
Then, the operator S can be represented as: 

5 5 1 

Sxy — Yq — ' with — — — S uv = - (S xu 5 yv + 5 xv 5 yu ) . (5-7) 

0O X y 0O X y Z 

The commutation relation (5.3) suggests the use of a Fock-space language, where 
and S are interpreted as creation and annihilation operators for dipoles in the projectile. 
More precisely, they create and, respectively, annihilate Pomerons, where a 'Pomeron' 
refers to the S'-matrix for an external dipole. Then, the Hamiltonian (5.4) can be viewed 
as an effective theory for Pomerons [11]. The first line of Eq. (5.4) describes the BFKL 
evolution of a Pomeron together with its dissociation : one Pomeron can split into two via 
the 'triple-Pomeron splitting vertex' (4.14), which is essentially the dipole kernel (2.3). 
The second line in Eq. (5.4) is similarly interpreted as the recombination of two Pomerons 
into one, via the 'triple-Pomeron merging vertex' (4.23) [33,34]. It is easily verified that 
the evolution equations for Pomeron correlations generated in this way, namely 

^ = (HS^) , (5.8) 

(where = SiS 2 ■ ■ ■ S K ) reproduce indeed the expected equations for the scattering am- 
plitudes, cf. Eqs. (2.1)-(2.2). Specifically, the Pomeron splitting piece of the Hamiltonian 
(the first line of Eq. (5.4)) generates the BFKL terms and the unitarity corrections, as 
originally noticed by Janik [35], whereas the Pomeron merging piece is responsible for the 
fluctuation terms [4,9,10]. 

But, clearly, this effective theory cannot be given a probabilistic interpretation — as a 
'reaction-diffusion process with Pomerons' — , for the same reason as for its dual counter- 
part discussed in Sects. 3.2 and 3.3: namely, the sign problem of the triple-Pomeron vertex 
for merging, Eq. (4.23). One may hope to circumvent this problem by reformulating the 
stochastic dynamics directly in terms of (projectile) dipoles, rather than the associated 
'Pomerons'. However, the subsequent analysis will show that this is not the case. 



5.2 Reaction- diffusion process in the projectile 



In this section, we shall study what should be the consequences — in terms of evolution 
equations for dipole scattering amplitudes — of a dipolar reaction-diffusion process taking 
place at the level of the projectile. That is, we shall assume that the projectile is a 
collection of dipoles which evolve via splitting and recombination with generic vertices, 
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and we shall use the evolution equations for dipole densities established in Sect. 4 to deduce 
equations for the scattering amplitudes between these dipoles and a generic target. 

In doing so, we shall follow a strategy pioneered by Kovchegov [27] in his construction 
of the Balitsky-Kovchegov equation [5,27], which is interesting in that it allows one to 
free oneself from the two-gluon exchange approximation: the obtained equations for the 
scattering amplitudes include multiple scattering to all orders. This method has been 
subsequently used by Levin and Lublinsky to give a rapid derivation of the large-iV c 
version of the Balitsky-JIMWLK hierarchy [28, 46]. In these early applications, the dipole 
evolution was restricted to splitting. More recently, Levin and Lublinsky have attempted 
to also include dipole recombination [12], but their final equations are plagued with some 
errors of sign that we shall correct in what follows. 

Before we proceed with our analysis, let us stress once again that the evolution to be 
described below does not represent the actual, physical, evolution of the projectile in 
QCD at large- N c (not even effectively !), but is only a fictitious reaction-diffusion process, 
whose consequences will be compared to the actual QCD equations in Sect. 5.3. 

The generic Master equation describing the would-be evolution of a system of dipoles 
in the projectile has been already written down in Sects. 4.1 and 4.2. The only other 
ingredient that we need for the present purposes is a factorization formula relating the 
projectile-target scattering amplitude to the evolving distribution of dipoles in the pro- 
jectile. Following Kovchegov [27], we write this amplitude as 

oo - 

A(Y) = W dr K (-1) K+1 n^( Ul ...u N ; Y - Y ) T^( Ul ...u N ; Y ), (5.9) 

K=l J 

where Y is an arbitrary intermediate rapidity with < Y < Y, n^{Y — Y ) is the 
K-dipole density in a projectile which has been evolved over Y — Y units of rapidity, and 
T^(Yq) denotes, as usual, the amplitude for the scattering between a set of k dipoles and 
a target evolved over Yq units of rapidity. Note that, in writing Eq. (5.9), no restriction 
has been assumed on either the structure of the scattering amplitude for the individual 
dipoles, or on the number k of dipoles which scatter simultaneously. Thus, Eq. (5.9) 
includes multiple scattering to all orders, as announced. 

Requiring that the total amplitude be insensitive to the rapidity divider Y between the 
target and the projectile (i.e., to the choice of a Lorentz frame): 

then using the evolution equations for the dipole densities (cf. Sect. 4) and the fact that 
these densities are numerically independent, one eventually obtains an hierarchy for the 
scattering amplitudes of projectiles with a given number of dipoles [12,28,46]. 

As before, we shall separately consider the contributions due to dipole splitting and merg- 
ing. For the splitting part, the results of Sect. 4.1 together with the argument above imply 
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dT(x) 



dY 

u 

and furthermore 

dT ( - 2 \x 1 x 2 ) 



dY 



a(x\uv) \-T(x) + T(u) + T(v) - T (2) (uv) 



J a( Xl \uv) [-T (2) ( Xl x 2 ) + 2T (2) (vx 2 ) - T (3) (uvx 2 ) 

uv 

+ \xi <-> x 2 y 



(5.11) 



(5.12) 



Higher equations in this hierarchy can be similarly obtained [28] 



The physical interpretation of these equations becomes more transparent if one recalls that 
a projectile dipole is more properly identified with the corresponding ^-matrix S = 1 — T 
(which describes its survival probability), rather than with its scattering amplitude T. 
Replacing T by 1 - S in Eq. (5.11) yields 



dS(x) 



dY 



-I 

a •> 



a(x \uv 



-S(x)+S {2) (uv) 



(5.13) 



where we have used the fact that T^ 2 \uv) = (T(u)T(v)) together with the simplified 
notations S(x) = (S(x)) and S^^uv) = (S(u)S(v)) . The evolution described by this 
equation is then interpreted as follows : When increasing Y by dY, the original dipole 
x can either split into a pair of dipoles (uv) with a probability density a(x\uv)dY, or 
remain in its original state with a reduced probability 1 — XdY, where A = J uv a(x\uv) is 
the inclusive splitting rate. 

Now, let us assume that the dipoles in the projectile can also recombine, with some 
generic vertex (3. Then, the evolution of the scattering amplitudes acquires additional 
contributions, which can be inferred from Eqs. (4.20) and (4.21) for the dipole densities, 
together with Eqs. (5.9)-(5.10). One thus finds: 



0T^( Xl x 2 ) 



dY 



P(x 1 x 2 \v) \T( Xl ) + T(x 2 ) - T(v) - T^( Xl x 2 ) 



(5.14) 



Once again, the physical interpretation becomes more transparent after replacing T = 
1 — S, which yields 



X\X 2 , 



dY 



J (3(x lX2 \v) [-S^( Xl x 2 )+S(v) 



(5.15) 



Namely, the two incoming dipoles x\ and x 2 can either recombine with each other into 
a dipole v with probability density f3(x\x 2 \v)dY , or survive in the original two-dipole 
configuration with the reduced probability 1 — dY J v f3{xix 2 \v). 

Let us finally notice an error of sign in the corresponding calculation in Ref. [12]: the 
analog of Eq. (5.14) - - that is, the contribution of the recombination process to the 
evolution of the scattering amplitude — appears in Ref. [12] with a sign which is opposite 
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to ours 10 . In view of the limpid physical interpretation of this equation given above, it is 
quite clear that the sign appearing in Ref. [12] is in fact wrong. 

5.3 Comparison with the evolution equations in QCD 

In this section, we shall compare the equations for the scattering amplitudes derived 
previously from the Master equation to the actual evolution equations in QCD, and thus 
conclude that the latter cannot be reproduced by a (probabilistically meaningful) dipole 
reaction-diffusion process within the projectile. 

Of course, there is no a problem at all in so far as the splitting alone is concerned: With 
the splitting vertex a identified with the dipole vertex in Eq. (4.14), the above equations 
(5.11) and (5.12) are immediately recognized as the first two equations of the Balitsky- 
JIMWLK hierarchy at large- N C) in agreement with the conclusions in Refs. [27,28]. 

Rather, the problem comes again from the recombination effects, as encoded e.g. in 
Eq. (5.14). When comparing the terms linear in T in the r.h.s. of this equation to the 
fluctuation term in the corresponding QCD equation (2.2), one comes across the same 
difficulty as previously discussed (towards the end of Sect. 4.2) in the context of target 
evolution. Namely, in order to reproduce the correct fluctuation term, one should iden- 
tify the recombination vertex (3 with the triple-Pomeron vertex in Eq. (4.23). With this 
choice, the first two terms in Eq. (5.14) would vanish (since the vertex (4.23) is a total 
derivative w.r.t. its last argument), while the third term there, which is surviving, would 
reproduce the fluctuation term in Eq. (2.2). But as repeatedly stressed, the vertex (4.23) 
has no fixed sign and thus cannot be used within a physical reaction-diffusion process. 
The incongruous nature of this choice becomes even more striking when one realizes that 
the would-be surviving term oc (—T(v)) in the r.h.s. of Eq. (5.14) yields a negative contri- 
bution in any physical recombination process (where /3 is positive semi-definite), whereas 
the fluctuation term that we need to reproduce here is rather positive ! Thus, the fluctua- 
tion terms in the equations for the scattering amplitudes cannot be reproduced, not even 
effectively, by a dipole recombination process in the projectile. 

Incidentally, the vertex (4.23) which is needed to formally relate the prediction (5.14) of 
the recombination process to the actual fluctuation term in QCD has the opposite sign as 
compared to the corresponding vertex used in Ref. [12] (see Eq. (3.30) there). The reason 
why the correct sign has been nevertheless reported in the final equations in Ref. [12] is 
because this sign difference in the vertex has been compensated there by a sign error in 
the equations for the scattering amplitudes, as mentioned in the previous subsection. 

Although the overall sign in front of the triple-Pomeron vertex (4.23) is conceptually 
unimportant — it does not affect our conclusion about the lack of a probabilistic inter- 

10 See Eq. (2.19) of Ref. [12], and focus on the terms involving 7„_i (the analog of our T^ -1 )) 
in the r.h.s. of the equation for j n (our T^). The vertex T2->i which appears there corresponds 
to our (3. Clearly, the terms proportional to 7 ra _i in the r.h.s. of Eq. (2.19) have the opposite 
sign as compared to our Eq. (5.14), or the more general equation (A. 7) in the Appendix. 
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pretation — , it is nevertheless interesting to clarify the sign difference between our vertex 
(4.23) and that in Eq. (3.30) of Ref. [12]. To that aim, we shall re-examine the argument 
used to obtain this vertex in Ref. [12] and show that, when carefully implemented, this 
argument leads in fact to the sign displayed in our Eq. (4.23). Perhaps the main interest 
of the subsequent discussion is to demonstrate that, also in the context of Ref. [12], the 
'dipole recombination' vertex has been obtained through a formal re-interpretation (sim- 
ilar to our above identification between the actual fluctuation term in Eq. (2.2) and the 
terms linear in T in Eq. (5.14)), which does not presuppose in any way the existence of 
an actual, physical, dipole recombination process. 

For definiteness, let us consider the fluctuation terms in the equations for the scatter- 
ing amplitudes. (By exchanging the roles of the projectile and the target, the following 
argument can be easily translated to the unitarity corrections.) As we have seen, these 
terms have indeed a dipole interpretation at large N c , namely, they correspond to dipole 
splitting in the target. Loosely speaking, splitting in the target corresponds to merging in 
the projectile (as this amounts to viewing the same Feynman diagrams upside down !), 
which makes it tempting to try and interpret the fluctuation terms as the result of dipole 
merging in the projectile. This in turn amounts to identifying the actual QCD process de- 
scribed in Fig. 1 to the fictitious process illustrated in Fig. 2, and then use this matching 
condition to extract an effective 'dipole recombination' vertex. As we shall shortly see, 
this procedure yields indeed the vertex shown in Eq. (4.23), including the sign. 

Specifically, the diagram on the r.h.s. of Fig. 1 is computed as 

(X f r i r 

Vl = 27 J M ^ U ' V ' ^ r Q & 4 °( a!l > Vl \ u > z n [-<x 2 sA)(x2, V2\z, v)\ + {1 -> 2}. (5.16) 



It is important to notice at this point that the elementary dipole-dipole scattering is 
i 2 a 2 Ao — — ct 2 A$ < 0, where the two factors of i come from the two propagators of the 
exchanged gluons. By assumption, the diagram in the r.h.s. of Fig. 2 is equal to 

V 2 = -a 2 J (3(x!, yi, x 2 , y 2 \x, y) A (x, y\u, v). (5.17) 

xy 



Setting Vi = V 2 , then acting with V^V 2 , where u ^ v, on both sides of the equation 
and using V^.ln(a; 2 ) = 4tt5(x), we finally obtain j3 = flsp, as anticipated. Note that 
the difference in sign as compared to the corresponding calculation in Ref. [12] comes 
from the minus sign in the structure of the dipole-dipole elementary amplitude (this sign 
has been overlooked in Ref. [12]). But even though it generates, by construction, the 
correct fluctuation terms, the 'recombination vertex' of Eq. (4.23) cannot be used within 
a meaningful stochastic process, because of the sign problem explained before. 

Ultimately, the lack of a probabilistic interpretation in terms of dipole recombination 
reflects the profound physical difference between the actual QCD mechanism for gluon 
saturation — the coherent suppression of gluon emission by strong color fields [31,44] - 
and the 'mechanical' 2 — > 1 dipole merging. 
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Fig. 1. Dipole splitting in the target 




Fig. 2. Dipoles (effectively) "recombining" in the projectile 
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A Evolution equations for a generic reaction— diffusion process 

In this section we present, for completeness, the general evolution equations for a generic 
reaction-diffusion process with non-local vertices for splitting and merging. Particular 
cases of these equations have already appeared in the main body of the paper. 

We consider the reaction-diffusion process A ^ A + A in which one particle at position x 
can split into two particles at positions x± and x 2 at a rate a(x\xix 2 ), and two particles 
at positions X\ and x 2 can merge into one particle at position x at a rate (5{x\X 2 \x). 
We denote by Pn{xi ■ ■ -Xn) = Pn(X) the probability density to observe a system of N 
particles at positions x±, . . . ,xn at a given time t (the dependence of Pn upon t is kept 
implicit). Due to merging and splitting, the probabilities evolve according to the following 
hierarchy of equations, which is usually referred to as the Master equation : 
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<9i-P/v(X)=^ / a(z\xiXj)P N -i(xi . . . fit . . . . . . x N z) - E / a(x i \z 1 z 2 )P N (X) 

ijtj Jz i Jzi ' Z2 

+E / P(ziz 2 \xi)P N+ i(xi . . . ^. . . . x N ziz 2 ) - E / P(xiXj\z)P N (X).(A.l) 

i Jz ^ i^j Jz 

In writing this equation, we have also assumed the symmetry properties a(x\xix 2 ) = 
a{x\x 2 X\) and fl(xix 2 \x) = j3{x 2 Xi\x) 

Let us now derive the equation for the evolution of an observable O with an average value 
defined by 

00 „ 

(C) = E / dF N P N (x 1 . . . x n )O n (xi . . . x N ), 

N=l 

where in the integral over the phase-space one needs to count only once the configurations 
differing only by the order of their elements (cf. Eq. (4.5)). After some straightforward 
algebra, we obtain the following evolution equation for the observable 



00 1 r N 

d t(°) = E 7777 M E / «M*i*2) \0 N+1 (x 1 ..4 i ...x N z 1 z 2 )-0 N (X) 

N=l Xi U=l Jz l> z 2 



+ E / P(Wj\z) [ON-ifa ..x N z) - O n (X)] \ ,(A.2) 

i+3 Jz J 



which encompasses Eqs. (4.8) and (4.19) of the main text. 

In particular, let us derive the evolution equations for the average A;-particle densities. 
The average particle number density (n(y)) is defined as in Eq. (4.9), and involves the 
following density operator 

N 

n N (y\xi . . .x N ) = E)<W ( A - 3 ) 
i=i 

Similarly, the /c-particle density operator is defined as 



n 



(y 1 ...y k \x 1 ...x N ) = <Wi 1 --- ( Wv 



Inserting the definition (A. 3) into the evolution equation (A. 2) gives 

d t (n(y))= / a{z x \yz 2 ) {n{z^)) + a{z 2 \z x y) (n(z 2 )) - a{y\z x z 2 ) (n(y)) (A.4) 



Zl,Z2 



+P(z lZ2 \y) {n (2 \ Zl z 2 )) -P(z iy \z 2 ) (n^(z u y)) - P(yz 2 \ Zl ) (n^(yz 2 ) 

In a similar way, but after some lengthier manipulations, one ends up with the following 
evolution equation for the general /c-particle density: 
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d t (n^( yi ...vy k )) (A.5) 
= / ^2a(z 1 \yz 2 ) (n (k \ yi . ..y i - 1 z 1 y i+1 . . .y k )\ - /3(z 1 y\z 2 ) (n {k+1 \ yi . ..y^)) 

+ a(z 2 \z 1 y) (n (k \y l . ..yi-^yt+i ■ ■ .yk)) - (3(yz 2 \z l ) (n {k+l \ yi . ..y k z 2 )) 
- a(y\z lZ2 ) (n {k \ yi ...y k )) + P(ziz 2 \y) (n {k+1) ( yi . . . f. . . . y k z x z 2 )) 
+ J^aiAviVj) (n {k ~ 1) (yi--^ i --^ r --y k z)) - P(yiyj\z) (n {k) ( yi . . .y k )) . 

Finally, we use the density operators to construct a scattering amplitude which includes 
multiple scattering to all orders. This is defined as [27] 

= E rr / (-) k+1 T^( Xl ...x k ; t ) n< fc >(*i ...x k ;t- t ) (A.6) 

k K \ Jxi 

with being the scattering amplitude for a set of k particles. In this expression to 
represents an arbitrary intermediate time, but the physical quantity A(t) is of course in- 
dependent of t : dA(t)/dt = 0. From this condition together with the evolution equation 
for the average densities derived before, one can finally deduce a set of evolution equations 
for the scattering amplitudes. The general equation in this hierarchy reads: 



d t T {k \x l ...x k ; t)=jr f a(x t \z lZ2 ) \t^\ Xi x kZl ; t) + T {k \ Xl . . . . . . x k z 2 ; 

i=1 Jz-L,Z2 

-T { - k \x l ...x k ;t)- T { - k+1 \x l . . . £ . . . y kZl z 2 ; t) 
+E / [T^Hx, ...f i ...x k ;t) + T^ k ~ 1 \x 1 .../ r ...r,: t) 

-T^" 1 ) (./•,...,/• y k z; t) - ( Xl ...x k ; t) 

(A 



i+3 
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